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Abstract 

A calculation of vector and tensor polarization observables of a deuteron in an optical 
potential framework for elastic proton-deuteron scattering at the incident deuteron 
energy Ed = 880 MeV is presented. The main aim is to point out the role of 
^ ■ various effects in explaining the wide-range angle distribution of these observables. 
Q ■ Under the investigation are a double-scattering mechanism in the optical potential, 
the influence of various iViV-partial waves and as well a comparison of calculations 
! with two different deuteron wave-functions derived from the Bonn-CD iViV-potential 
model and dressed bag model of Moscow- Tuebingen group. 
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(N 1 Introduction 

The reaction of elastic nucleon-deuteron scattering are considered both by experimen- 
talists and theoreticians as one of the clue problems in few-nucleon physics. For three 
decades it has been served as a hope to obtain more information about intermediate- and 
short-range NN interaction and as a probe of the deuteron structure at small distances 
(for a review, see Ref. []]). During the last decade it is also studied with the purpose 
of testing of various 3iV-forces and particularly their spin-dependance [21 E] • It is also of 
| interest as a basic reaction to establish a polarimetry for vector-tensor mixed polarized 
^ ' beams. 

Now, there is a set of high quality nucleon-nucleon potentials that describe two-nucleon 
elastic scattering data perfectly below the pion production threshold. But they underesti- 
mate the three nucleon bound energy. The incorporation of a 3N force in the calculations 
improves the situation but some unsolved questions remained, especially it concerns its 
spin-dependant properties. The study of polarization observables in elastic pd scattering, 
as expected, can provide an additional information about the properties of 3N force. To 
describe pd elastic scattering below the pion production threshold different techniques have 
been applied. Among them the momentum space Faddeev equations can now be solved 
with high accuracy for the most modern two- and three-nucleon forces. The 3N forces in 
such a calculations are of Fujita-Miyazawa [U [5] or Tucson-Melbourne [6] type. Also, in 
Ref. [7j a coupled- channel model with excitation of the stable A-isobar was successfully 
applied to pd scattering employing Chebyshev expansion of two-body amplitudes. 

But calculations at incident proton energies greater than 200 MeV in lab encounter 
with some nontrivial difficulties. The first problem is that up to now there is no reli- 
able quantitative model for iVW-interaction above the inelastic threshold. All existing 
models that pretend to description of two-nucleon scattering up to 1 GeV do it only in 
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semi-quantitative manner [8]. The second problem is that it is no longer conventional 
three-body problem above the pion production threshold. The nucleon isobar and meson 
degrees of freedom start to play a significant role. And last but not least, there are purely 
computational difficulties in performing relativistic three-body calculations at higher en- 
ergies [;9|. Hence one must resort to some kind of simplifications and approximations in 
the study of this reaction. 

In this work, the optical potential framework to calculate polarization observables 
of a deuteron is used. The incident deuteron energy is taken to be 880 MeV in lab. 
This energy is just below the region where the influence of A-isobar excitation becomes 
significant masking a large momentum component of the deuteron wave function [TO] . 
Moreover, this energy lies in the interval of energies proposed for the measurements of 
dp elastic scattering analyzing powers at internal target station of the Nuclotron at JINR 

P-D- 

2 Theoretical framework 

The basic equations for three-body problem in quantum physics are the Faddeev equations 
which can be written in operator form: 

U = PG 1 + PTGqU, (1) 

where U = ^ Mm (q',q) is an amplitude of elastic dp-scattering, fi d , fi N are spin 
quantum numbers, T - NN scattering matrix, Go = (E — H Q + ie)' 1 is a free propagator 
of 3iV system and P = P12-P23 + -P13-P23 stands for a permutation operator that takes into 
account the property of identity of three nucleons. Here and further in this work we have 
dropped out all 3N force contributions. 

One can rewrite (CQ) in the form that is more appropriate for calculations at higher 
energies where an impulse approximation is justified: 

U = V opt + V opt G d U. (2) 

This is the optical potential framework and V opt is a nucleon-deuteron optical potential: 

V apt = PGo 1 + PT G V apb . (3) 

Here T c is a NN T-matrix with subtracted deuteron pole term and Gd is a deuteron 
contribution in the spectral decomposition of the two-body Hamiltonian. 

Although, the Eq. ([2]) seems quite simple and represents a common two-body scat- 
tering equation, the derivation of the optical potential by solving Eq. (jSJ) is as difficult 
as solving Faddeev equations themselves. However, as was mentioned in Introduction, at 
energies higher than the pion production threshold, new degrees of freedom, other than 
nucleons, come into the game. But this new physics is not taken into account in (j3j). 
Also, these are the energies at which the impulse approximation is a good starting point, 
so one can make a series expansion of V opt in ([3]), keeping a possibility to include in V opt 
some additional terms: 

K P t = PG 1 + PT C P + PT C G PT C P + ... 
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Figure 1: Pictorial representation of pd optical potential up to the second order in the 
series expansion of Eq. (j3J) . 



Pictorially it can be described as the sum of nine graphs (up to the second order in series 
expansion), see Fig. [TJ The first diagram is a usual one-nucleon exchange mechanism, the 
next three diagrams are so called triangle diagrams and represent the single-scattering 
process and the last five diagrams take into account the inelastic break up channel. 
The one-nucleon exchange contribution is written in a usual way: 
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Here, we omit a complicated notation related to a summation over spin and orbital quan- 
tum numbers. ^23 is a wave function of the deuteron composed of nucleons 2 and 3. 

To evaluate a single scattering diagram, one must implement an integration over the 
internal momentum of nucleons in the deuteron. To do this, a knowledge about the fully 
off-shell behavior of NN T-matrix is required: 
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here k = q — q' is a transferred momentum. 

However, the absence of high quality NN interaction model above the inelastic thresh- 
old at the present moment forces anyone to make some approximate evaluation of the 
integral, proceeding from the assumption of either an off-shell behavior of somehow pa- 
rameterized T-matrix or taking only its on-shell value. In this work, the T-matrix is put 
on-shell in the way that minimizes the off-shell corrections. This is an optimal impulse 
approximation [121 H3] : 



T c (q',p-jq' +\q\ q,P~jq +\q') — > T f q', -|ep +i<f; q, -|g +^q'; E eff 
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Then T-matrix is evaluated in the so called NN Breit frame, where an effective scattering 
energy is 
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q 2 + Mi N +\l^(5-3cos6) + Mj [ j -jcos 2 ^. 

For example, at deuteron laboratory energy Tj a b = 880 MeV and scattering angle 9 = 10° 
the effective scattering energy of two nucleons is E e g = 443 MeV, whereas at the angle 
9 = 170° - E e Q = 796 MeV. The numerical values of the T-matrix are obtained from 
partial wave analysis SAID [13] that provides the NN phase shifts up to 3 GeV and the 
total angular momentum J NN = 7. 

This permits the T-matrix to be shifted outside the integral and the remaining inte- 
gration produces the deuteron form factor: 

i(23)(/|T c |z) = T c (V, -\q> +\q- q, -\q+\q>; E e// ) S(k). 

For the double scattering term one has two integrals - one over the internal momentum 
in the deuteron and another over the intermediate momentum of the scattered nucleon: 



1(23) < 



x ^ 23 (t=E -Zq* -\q +1/) T 2 (q'J -\q' +±g 7 ; g 7 ',p' -fg 7 +±g 7 ) G 

x Tx (g 7 p -fg 7 +\q; q,p -fg +±g 7 ') ^ 23 (^+|g +|g 7 -\i 
Taking into account in Go only its pole part and employing the impulse approximation 
T2 (q',p' -§<? +\q ! '\ q 7 ',p' -fg 7 +\q) GqT, (g",p -fg 7 ' +\q; q,p -\q +±g 7 ) — > 

- - ?)T 2 (g', -fg 7 +±g 7 '; g 7/ , -fg 7 +±g 7 ; E eff2 ) Tx (g 7 , -fg 7 +±g; g, -fg-+±g 7 ; £ e //i) , 

finally one gets 

i(23)(/|T 2 Go7 1 i|i) = -IvrgMjv J ^^T 2 (q> , -\q> +±g 7 ; g 7 , -fg 7 +±g 7 ; E eff2 ) 

Tx (g 7/ , -fg 7 ' +|g1 g, -|g+^ 7/ ; E eff i) s (j> -\{q' + g)J • 

The scattering equation (T5]) is solved in helicity basis employing f^-matrix approxima- 
tion, i.e. only the pole part of the two-body propagator Gd is remained thus all terms in 
the equation contain only on-shell information about the optical potential and scattering 
amplitude: 

(A^l?/- 7 ^, g)|A d A^> = (A'^lV^tC?', 9)IA d A^>-^^^<A' d A / Ar IV;^t(?^ ^^(g, g)|A d A,v>+ 
l V J -qT^i (< WlK J pt (g', q"W\q\ q)\X d X N )A(q", qj-^'MV^q ', q)U J (q, q)\X d X N )A(q, g)) , 
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where the kinematical factor is 
A(q", q) = q 



m (E q + E q „){{El + El„)/2 - q* - g" 2 - M% - Mj] 



El + E\„ 

Then for the scattering amplitude in spin space one has (the incident particle is going 
along the z-axis): 

x d J x , N _ x>N+ ^ d (6)(\' d \' N \U J (q,q)\ -Hdfi N ). 
It was found that at = 880 MeV a convergent result is obtained at J = 37/2. 



3 Results 

The calculation of the polarization observables is performed with to different deuteron 
wave functions. The first one is a wave function derived in the meson-exchange Bonn- 
CD model [15]. The general trait of wave functions of this kind, derived from the most 
modern NN potentials, is their depletion at small distances. The other possible choice is 
the wave function with a nodal behavior. This node corresponds to a so-called forbidden 
state in NN system as a consequence of the six-quark dynamics and the fact that the 
mostly symmetric six-quark state |s 6 ) has a small NN component [16J. As one of the 
representatives of such a wave functions with a nodal behavior can serve the wave function 
of dressed bag (DB) model of Moscow- Tuebingen group [17] . In Fig. [2] a comparison of 
two calculations with these different wave function is presented. As one can see, the 
remarkable difference in the polarization observables occurs at large scattering angles at 
which the high momentum component of the wave functions is probed. It should be 
noted however, that even if we can discriminate this two choices of the wave functions by 
experimental data, to obtain an ultimate conclusion, which of the models are preferable, 
one has to include in calculations a corresponding 3N force. And the 3N forces are 
different in these two concepts. 

In Fig. [3] the influence of rescattering and break up channel on the polarization 
observables at the deuteron energy 880 MeV is shown. It is evident that taking into 
account the break up process or double-scattering in the optical potential is inevitable 
to obtain quantitative predictions at these energies. The rescattering terms which are 
imbedded in the pd scattering equation is not so pronounced and much less important, 
however not negligible. 

Finally, the contribution from different partial waves in iViV T-matrix is examined, see 
Fig. HI The full calculation takes into account all NN partial waves that parameterized 
in SAID partial wave analysis [H], i.e. Jnn = 7. The other two calculations are done 
with total angular momentum up to Jnn = 5 - dashed curve and J^n = 3 - dotted curve. 
Although at some scattering angles the differences between the observables are very small, 
the convergence is only partial. But it seems that Jnn — 7 is sufficient angular momentum 
at this energy range. 
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4 Summary and Conclusions 

A calculation of the deuteron polarization observables Ay, A yy , A xx and A xz through an 
optical potential model for elastic proton-deuteron scattering at the incident deuteron 
energy = 880 MeV was presented. Under investigation were the double-scattering 
mechanism in the optical potential, the effect from precise treatment of two-body unitarity 
in proton-deuteron Lippman-Schwinger equation, the influence of various iViV-partial 
waves and as well a comparison of calculations with two different deuteron wave-functions 
derived from the Bonn-CD iViV-potential model (15] and DB model of Moscow- Tuebingen 
group [IT]. For the NN input, the model independent approach in which nucleon-nucleon 
T-matrix is taken to be on-shell and evaluated in the so-called Breit frame, was used. The 
calculation was carried out in the helicity partial wave decomposition and total angular 
momenta up to J = 37/2 in proton-deuteron system seem to be sufficient to obtain 
a convergent result. As for iViV-partial waves, we took all SAID partial waves up to 
Jnn — 7. The large effect from double-scattering mechanism at 880 MeV as well as 
the large dependency on the short-range behavior of the deuteron wave function were 
observed. 

This work is partly supported by the Russian Foundation for Basic Research, grant 
04-02-17107a. 
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Figure 2: Deuteron vector and tensor polarization observables at the energy Ed = 880 
MeV in lab. The solid and dashed curves are calculations with the deuteron wave function 
derived in the Bonn-CD model and DB model of NN interaction respectively. 
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Figure 3: The effect of a break up channel in the optical potential - dashed curve, and 
when all rescattering mechanisms are omitted - dotted curve. The solid curve is a full 
calculation with Bonn-CD deuteron wave function. 
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Figure 4: The results of calculations that take into account different amount of partial 
waves in NN T- matrix. The solid curve is a result when all J NN = 7 SAID partial waves 
are included in calculation, whereas dashed and dotted curves are the calculations with 
only Jnn = 5 and Jnn = 3 respectively. 
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